Sys.setenv(TZ="UTC")
t0 <- as.POSIXct("2023-07-04")
call <- knitr::opts_current$get("params.src")
dpath <- if (is.null(call)) "." else "resources"
'getTrack' <- function(device,vehicle=FALSE,ltraj=FALSE,inside=NULL) {
   if (missing(device))
      return(NULL)
   if (!ursa::is_spatial(device)) {
      list1 <- dir(path=file.path(dpath,device)
                  ,pattern=".*20[x2]3.+\\.gpx$",full.names=TRUE)
      res <- lapply(list1,\(fname) {
         a <- sf::st_read(fname,layer="track_points",quiet=TRUE)[,1:5]
         a <- a[!is.na(a$time),]
         if (!ursa::spatial_count(a))
            return(NULL)
         taf <- table(track_fid=a$track_fid)
         tas <- table(track_seg_id=a$track_seg_id)
         if (length(tas)>length(taf)) {
            b <- by(a,a$track_seg_id,\(x) x)
         }
         else {
            b <- by(a,a$track_fid,\(x) x)
         }
         a <- lapply(b,\(x) {
            rt <- range(x$time)
            y <- ursa::trackline(x)
            len <- ursa::spatial_length(y)
            dur <- as.numeric(diff(range(x$time)),units="hours")
            speed <- len*1e-3/dur
            sa <- data.frame(n=ursa::spatial_count(x),length=len,dur=dur,speed=speed
                            ,from=rt[1],to=rt[2])
            cond1 <- sa$dur>10/60 & sa$n>10
            sa <- sa[cond1,]
            if (!cond1)
               return(NULL)
            x$device <- tolower(device)
            x$segment <- digest::digest(x,"crc32")
            tail(x[,c("time","ele","device","segment")],-1)
         }) |> do.call(ursa::spatial_bind,args=_)
         a
      }) |> ursa::spatial_bind()
      res <- res[order(res$time),]
      if (vehicle) {
         if (length(ind <- grep("^(cruise|road)",res$survey,invert=TRUE))) {
            res <- res[-ind,]
         }
      }
      if (!is.null(inside)) {
         cond <- vector("list",nrow(inside))
         for (i in seq(nrow(inside))) {
            if (!inside$fill[i])
               cond[[i]] <- integer()
            else {
               j <- which(res$time>inside$begin[i] & res$time<inside$end[i])
               cond[[i]] <- j
            }
         }
         cond <- sort(unlist(cond))
         if (!length(cond))
            return(NULL)
         res <- res[cond,]
      }
      res$id <- "merge"
   }
   else {
      res <- device
      ltraj <- TRUE
   }
   if (!ltraj)
      return(res)
   if (length(ind <- which(duplicated(res$time)))) {
      res <- res[-ind,]
   }
   res2 <- adehabitatLT::as.ltraj(sf::st_coordinates(res)
                                 ,date=res$time
                                 ,id=res$id
                                 ,proj4string=sp::CRS(sf::st_crs(res)$proj4string)
                                 )
   opW <- options(warn=-1)
   res3 <- adehabitatLT::cutltraj(res2,"dt>c(days=1)*c(hour=1)*c(min=12)*c(sec=60)"
                                 ,nextr=TRUE)
   options(opW)
   res3
}
'inside' <- function(device,dt=as.difftime(12,units="mins"),dl=50) {
   if (is.null(device))
      return(NULL)
   ind1 <- which(diff(device$time)>dt)
   begin <- c(t0,device$time[ind1],max(device$time)+1)
   end <- c(min(device$time)-1,device$time[ind1+1],as.POSIXct(Sys.Date()+1))
   ret <- data.frame(from=c(NA,ind1,NA),to=c(NA,ind1+1,NA)
                    ,device=head(device$device,1)
                    ,begin=begin
                    ,end=end
                    ,length=NA
                    ,fill=TRUE
                    )
   for (i in seq(nrow(ret)) |> sample()) {
      ind <- c(ret$from[i],ret$to[i])
      if (anyNA(ind))
         next
      ret$length[i] <- ursa::trackline(device[ind,]) |> ursa::spatial_length() |> round()
   }
   ret$dt <- as.numeric(ret$end-ret$begin,units="hours") |> round(1)
   ret$fill[ret$length<dl] <- FALSE
   ret
}
'suppl' <- function(main,suppl,vehicle=FALSE) {
   more <- getTrack(suppl,inside=inside(main),vehicle=vehicle,ltraj=!TRUE)
   main <- ursa::spatial_bind(main,more)
   main <- main[order(main$time),]
   main
}
priority <- list(oregon=c("Oregon","Zhenya"))
pname <- names(priority)
for (i in seq_along(priority)) {
   message(pname[i],":")
   vehicle <- pname[i] %in% "general"
   main <- getTrack()
   for (device in priority[[i]]) {
      message("   ",device,":")
      main <- suppl(main,device,vehicle=vehicle)
   }
   print(getTrack(main,ltraj=TRUE))
   main <- main[,c("id","time","ele","device")]
   fname <- file.path(dpath,paste0("track-",pname[i],".geojson"))
   ursa::spatial_write(main,fname)
   tr <- ursa::trackline(main,by=main$id,connect="united")
   dst <- file.path(dirname(fname),gsub("^track","trackline",basename(fname)))
   ursa::spatial_write(tr,dst)
}

require(ursa)
(trkpt <- spatial_read("input/track-general.geojson"))